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We carry out a study of intermediate-mass (between 1 and 2.5 GeV) dilepton spectra from 
hadronic interactions in heavy-ion collisions. The processes considered are nn — > //, np — > 11, 
na\ —* 11, ncj — > 11, KK — > 11, and KK* + c.c —* 11. The elementary cross sections for those 
are obtained from chiral Lagrangians involving pseudoscalar, vector, and axial-vector mesons. The 
respective electromagnetic form factors are determined by fitting to experimental data for the reverse 
ON . processes of e + e — > hadrons. Based on this input we calculate cross sections and thermal dilepton 

ON ■ emission rates and compare our results with those from other approaches. Finally we use these 

elementary cross sections with a relativistic transport model and calculate dilepton spectra in S+W 
collisions at SPS energies. The comparison of our results with experimental data from the HELIOS-3 
collaboration indicates the importance of the secondary hadronic contributions to the intermediate- 
mass dilepton spectra. 
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The experimental measurement and theoretical investigation of dilepton production constitute one of the most 
active and exciting fields in the physics of relativistic nuclear collisions W . Because of their relatively weak final-state 
interactions with the hadronic environment, dileptons, as well as photons, are considered ideal probes of the early 
stage of heavy-ion collisions, where quark-gluon-plasma (QGP) formation is expected (|||. Because of an additional 
variable, the invariant mass M^, dileptons have the advantage of a better signal to background ratio than real photons 
They of course also prove superior in processes involving two-body annihilations. 
Dilepton mass spectra produced in heavy ion collisions can basically be divided into three regions. The low-mass 
£^ . region below (~ 1 GeV) is dominated by hadronic interactions and hadronic decays. In the intermediate-mass 
region between and about 2.5 GeV, the contribution from the thermalized QGP might be seen ^-Q). In the high- 
• rH , mass region at and above mj/q, the major effort in heavy ion experiments has been the detection and understanding 
X ' of J/* suppression. 

So far, the experimental measurement of dilepton spectra in ultrarelativistic heavy-ion collisions has mainly been 
carried out at the CERN SPS by three collaborations: the CERES collaboration is dedicated to dielectron measure- 
ments in the low-mass region ||||, the HELIOS-3 JlCj collaboration has measured dimuon spectra from threshold up 
to the J/\l/ region, and the NA38/NA50 Jll[] collaboration measures dimuon spectra in the intermediate- and high- 
mass re gions, emphasizing J/^S suppression (for a summary of low- and intermediate-mass dilepton measurements see 
Refs. P Jil^ , |l"3|l ) . In addition, dilepton spectra in heavy-ion collisions at energies of a few GeV/nucleon were measured 
by the DLS collaboration pi] . In the near future, dilepton spectra will be measured by the PHENIX collaboration 
|| at RHIC, and by the HADES collaboration at the GSI |§. 

Recent observation of the enhancement of low-mass dileptons in central heavy-ion collisions by the CERES |^|,^| and 
the HELIOS-3 |[(| collaborations has generated a great deal of theoretical activity. Different models have been used 
to interpret these data. The results from many groups with standard scenarios (i.e., using vacuum meson properties) 
are in remarkable agreement with each other, but in significant disagreement with the data: the experimental spectra 
in the mass region from 0.3-0.6 GeV are substantially underestimated 0,18| (see also Ref. j9|). This has led to 
the suggestion of various medium effects that might be responsible for the observed enhancement. In particular, the 
dropping vector meson mass scenario [ fi7|,|l9[|20f| is found to provide a unified description of both the CERES and 
HELIOS-3 data. However, see also Ref. f2j. 

In the high-mass region around mj/y, the J/ty suppression has been a subject of great interest, since it was first 
proposed as a signal of the deconfinement phase transition 122]. Various investigations show that up to central S+Au 
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collisions, the normal pre-resonance absorption in nuclear matter can account for the observed J/^S suppression p3[ . 
However, recent data from the NA50 collaboration for central Pb+Pb collisions show an additional strong 'anomalous' 
suppression which might indicate the onset of the color deconfinement 24 1. 

Other interesting experimental data that have not yet received much theoretical attention are dilepton spectra in the 
intermediate-mass region from about 1 GeV to about 2.5 GeV. Both the HELIOS-3 and NA38/NA50 collaborations 
have observed significant enhancement of dilepton yield in this mass region in central S+W and S+U collisions as 
compared to that in proton- induced reactions (normalized to the charged-particle multiplicity) pjJilfl. Preliminary 
data from the NA50 collaboration also show significant enhancement in central Pb+Pb collisions |11| (see also Ref. 

§)• 
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FIG. 1. Comparison of background with experimental data in p+W and S+W collisions. The data is from Ref. jioj 



For dilepton spectra with mass above 1 GeV, the contributions from charm meson decay and the initial Drell-Yan 
processes begin to play a role. These hard processes scale almost linearly with the participant nucleon number, 
and can therefore be extrapolated from proton-proton and proton-nucleus collisions. Such a study has recently been 
carried out by Braun-Munzinger et al [M. The results for p+W and central S+W collisions corresponding to the 
HELIOS-3 acceptance are shown in Fig. |l|, and are taken from Ref. These, together with the dileptons from the 
decay of primary vector mesons, are collectively termed 'background' in this work. It is seen that this background 
describes very well the dimuon spectra in p+W reactions, shown in the figure by solid circles. 

However, as can be from the figure, the sum of these background sources grossly underestimates the dimuon yield 
in central S+W collisions, shown in the figure by open circles. Since the dimuon spectra are normalized by the 
measured charged particle multiplicity, this underestimation indicates additional sources to dilepton production in 
heavy-ion collisions. There are at least three possible sources for this enhancement: the additional production of 
charmed mesons and/or Drell-Yan pairs, a QGP formed in the collisions, and secondary hadronic interactions. While 
all these possibilities are of interest, and may actually coexist, in this work we concentrate on the contributions 
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from the secondary hadronic interactions, which we believe need to be quantitatively assessed. In this work we limit 
ourselves to meson interactions. However we will also comment on the role of baryons later. For dilepton spectra 
at low invariant masses, it is well known that the tttt annihilation plays an extremely important role in heavy-ion 
collisions. It is also expected that other secondary processes will play a role in the dilepton spectra in the intermediate 
mass region, and we attempt to demonstrate this in a realistic calculation that compares with experimental data. 
This will be done with the relativistic transport model used in Ref. 0,^6) for the study of low-mass dilepton and 
photon production. The main motivation for such a study is to understand the origin of the observed enhancement 
in the intermediate-mass region, and to see whether the data calls for the formation of a QGP phase. 

In Section 2 we discuss the elementary dilepton production cross sections that are needed as inputs in the transport 
model. We shall emphasize the constraints imposed on these cross sections by the experimental data for the reverse 
process of e + e _ — > hadrons. In Section 3, we discuss the thermal dilepton emission rates from the hadronic interac- 
tions, and compare our results with those from other approaches. The elementary cross sections, as constrainted by 
the e + e~ annihilation data, are then used in the transport model to calculate dilepton spectra in heavy-ion collisions 
at SPS energies. The results will be presented in Section 4. The paper ends with a brief summary and outlook in 
Section 5. 



II. CROSS SECTIONS AND FORM FACTORS 



In the low-mass region from 2m w to m„o >ul , it has been shown that ir + Tr~ annihilation, which is characteristic 
of heavy-ion collisions, plays an important role in explaining the observed enhancement. Similarly, it is expected 
that other secondary hadronic processes should also play some role in the intermediate-mass region. Indeed, previous 
thermal rate calculations based on kinetic theory show that in the mass and temperature regions relevant for this study, 
the following processes (from the hadronic phase) are important: irir — > 11, np — ► 11, nai — > 11, ttlu — > 11, KK — » 11, 
and KK* + c.c — > U [|27|— 130|] . Am ong th ese processes, nai — > U has been found to be the most important one, mainly 
because of its large cross section |28| , |30|] (similar conclusions have been drawn for thermal photon production f3l"| , |32|| ) . 
Note here that we discuss only two-body reactions, as they should dominate the phase space region we are considering. 

An important input in the transport model calculation of dilepton spectra in heavy-ion collisions is the elementary 
dilepton production cross sections for the processes outlined above. In next three subsections, we will discuss the 
cross sections and form factors for pseudoscalar-pseudoscalar, pseudoscalar-vector, and pseudoscalar-axial vector 
meson annihilation. 



A. pseudoscalar-pseudoscalar annihilation 




FIG. 2. Feynman diagram for the annihilation of two pseudoscalar mesons into a lepton pair in the VMD model. 



In this class we have it — it and K — K annihilation. In the Vector Meson Dominance model (VMD), the Feynman 
diagram for this process is quite simple and is shown in Fig. @, where P represents a pseudoscalar meson, V the 



intermediate vector meson, 7* the virtual photon, and 11 the lepton pair. See Ref. 1 33 1 for a discussion of the possible 



uncertainties in VMD. We use a standard Lagrangian for the pseudoscalar-pseudoscalar-vector interaction, 

C V pp = gvppV^P 9 M P (1) 
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The obtained dilepton production cross section in pion-pion annihilation is well known |34| , |35| ] to be 

/ -i- ,tn 8na 2 k. 4m, . . 2m,. 

- «) = W'^ (M)I (1 - M^ ){1 + M^' (2) 

where k is the magnitude of the three-momentum of the pion in the center-of-mass frame, M is the mass of the lepton 
pair, and mi is the mass of the lepton. It is well known that the electromagnetic form factors l-F^M)! 2 play important 
role in this process, providing empirical support for VMD: the pion electromagnetic form factor is dominated by the 
p(770) meson. In addition, at large invariant masses, higher p-like resonances such as p(1450) were found to be 
important pG|. 

Using Eq7(0) and detailed-balance, we can get the cross section for 7r + 7r~ production in e + e~ annihilation, 

- ( e + e-- 7 r + ^) = ^^|F 7r (M)| 2 (3) 

This cross section has been measured with high accuracy in the mass region that is relevant to this study |3^ , |3^] . In 
Rcf. p(| a detailed analysis of the experimental data was carried in order to determine the pion electromagnetic form 
factor. Four p-like vector mesons were found to be present. The comparison of the form factor determined in Ref. |3(| 
with the experimental data from the OLYA collaboration |37| (circles), the CMD collaboration (squares), and 
the DM2 collaboration |Q (triangles) is shown in Fig. [| A direct comparison with the experimental cross section 
is shown in Fig. ^. The agreement between the model and the experimental data is excellent. This assures us that 
the elementary dilepton production cross section in pion-pion annihilation used in our transport model is well under 
control. 




FIG. 3. The pion electromagnetic form factor. The solid curve is based on the model of |3q| , while the dotted curve is based 
on the parameterization of one of us with Kapusta Q|. The experimental data are from the OLYA collaboration |^7| (circles), 
the CMD collaboration j37| (squares), and the DM2 collaboration (38) (triangles). 
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The pion electromagnetic form factor has been parameterized in different ways. For example, in Ref. |34|], one of 
us and Kapusta proposed the following form, 



in 



{M 2 -m' 2 ) 2 + (m r T r ) 2 ' 



(4) 



where m r — 0.775 MeV, m' r — 0.761 MeV, and T r = 0.118 MeV. This parameterization reproduces the Gounaris- 
Sakurai |3!| formula and describes very well the experimental form factor from the 2ir threshold to about 1 GeV in 
invariant mass, as shown in Fig. | by the dotted curve. It however underestimates the data at larger invariant masses, 
since it neglects the role of higher p-like resonances. 
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FIG. 4. The cross section for e + e — » n + n . The solid curve is based on the model of |3(J. The experimental data are from 
the OLYA collaboration (3t| (circles), the CMD collaboration (squares), and the DM2 collaboration H| (triangles). 



The dilepton production cross section in kaon-antikaon annihilation is very similar, 



(j{K + K-,K°K° -> 11) = 



8na 2 k 
3M 3 



\F K+tK0 (M)\ 2 (l 



4mf 
'M 2 
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2mf 
"M 2 



(•5) 



where k is the magnitude of the kaon momentum in the center-of-mass frame. The cross section for the reverse process 
of electron-positron annihilation can again be obtained from the detailed-balance relation, 



(6) 



In these equations, \F^+ | 2 and \Fk° | 2 are the electromagnetic form factors of charged and neutral kaons, respectively. 
These form factor are dominated by the phi meson, 0(1020). At higher masses, other vector mesons may become 
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important |36|]. The experimental form factor will be used in this work. It is shown in Fig. The solid and 
dotted curves are for charged and neutral kaons respectively, and the symbols are the experimental data for the 
charged kaon form factor from the CMD-2 collaboration [|o) (circles), the DM2 collaboration (squares), and the 
OLYA collaboration 0] (triangles). The neutral kaon form factor follows from general arguments |pq| . The direct 
comparison with the experimental charged kaon cross section is shown in Fig. ^. Again, very good agreement between 
the model and the data is seen. 
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FIG. 5. The kaon electromagnetic form factor. The solid and dotted curves are for the charged and neutral kaons, respectively 
[^6| . The symbols are the experimental data for the charged kaon from the CMD-2 collaboration ^(J (circles), the DM2 
collaboration (squares), and the OLYA collaboration jl2| (triangles). 
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FIG. 6. The cross section for e + e~ — > K + K~ . The solid curve is based on the model of |36). The experimental data are 
from the CMD-2 collaboration pCI (circles), the DM2 collaboration BlJ (squares), and the OLYA collaboration wA (triangles). 



B. pseudoscalar- vector meson annihilation 

In this class we include the following processes, 7rp — ► 11, KK* + c.c. — > 11, and wu> — » iZ. The first two processes 
effectively involve three pions, while the third one involves four pions. In transport models a process involving three 
or more pions in the initial state can only be described as a two step process with an intermediate resonance. The 
empirical success of transport approaches gives some credibility to this scenario. Note however that the treatment of 
quantum interference remains a possible issue in this framework |43[ . 

The interaction Lagrangian for this case is 

Cvvp = 9vppe^ af 3d^V v d a V f3 P . (7) 
The first two channels have been studied in Ref. [p9| . The cross section for 7rp annihilation is given by 

, I -, 27ra 2 fc 7r 4m? \ / 2m? \ , , 

where is the magnitude of the pion momentum in the center-of-mass frame. Note that the above cross section 
is evaluated in the narrow-width approximation for illustration purposes only. This simplification is not used in the 
actual transport calculation. The cross section for e + e _ — * np can again be obtained from detailed-balance. 

The electromagnetic form factor \Fi Tp (M)\ 2 can then be determined by analyzing the experimental data for e + e~ — > 
tt + tt~tt°. In Ref. pS[ |, three isoscalar vector mesons, 0(1020), u>(1420), and k>(1670) were found to be important in 
order to fit the experimental data from Refs. |44|j4^ ] (see also Ref. pq|), namely, 
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^w=i;(^) ^ tti: 4 ... , • (9) 



/ (my — M 2 ) — imyVy 



Here the summation runs over the three vector mesons listed above. While the coupling constants g^ and g^p can 
be determined from the measured widths, the coupling constants for other two mesons and the relative phases were 
determined by fitting to the experimental data of Refs. HJI^ j. In this work, we determine these coupling constants 
by fitting to the latest data from the DM2 collaboration ||47(] and the ND collaboration fig] . These parameters are 
only slightly different from those listed in Ref. p9| . The comparison of our fit with the experimental data is shown in 
Fig. @ 




M (GeV) 



FIG. 7. The cross section for e + e — » (tvp) — > ty + tt tv°. The solid curve is our fit to the experimental data from the ND 
collaboration Jtq] (circles) and the DM2 collaboration J47J (squares). 



Similarly, the cross section for KK* (or KK*) is given by 

«^^-^lW(i-£f)( 1 + ^), do, 

where kx is the magnitude of the kaon momentum in the center-of-mass frame. Since a K* eventually decays 
into a kaon and a pion, the electromagnetic form factor in the above equation can be determined by analyzing the 
experimental data for e + e~ — > K°K ± tt t , as was done in Ref. [^J. The form factor was found to be dominated by 
</>'(1680), 

TP (M\ - 9 ^' KR * nn 
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where the coupling constant 



9*1 



0.19 GeV 1 was determined by fitting the experimental data of the DM1 



collaboration J49J and the DM2 collaboration l5Cf] . The comparison with the experimental data is shown in Fig. 
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FIG. 8. The cross section for e e — > (KK* + c.c. 
experimental data are from the DM1 collaboration Wt 



-> i^iT^. The solid curve is based on the model of Ref. 
(circles) and the DM2 collaboration fed] (squares). 



The 



Finally, the cross section for lepton pair production in pion-omega annihilation is given by 

4na 2 k„ 



a{n { >uj -> 11) 



, / 4to?\ / 2mf\ 



9M 



(12) 



where fc^ is the magnitude of pion momentum in the center-of-mass frame. Experimentally, the electromagnetic form 
factor \Ft TUJ (M)\ 2 was studied in Ref. (48| by measuring e + e~ — > 7r°7r°7. The form factor was parameterized in terms 
of three isovector p-like vector mesons, p(770), p(1450), and p(1700) in Ref. p8[, 



F^{M)=Y J 



9v* 



9V 



(m 2 v 



M 2 



(13) 



Here the summation runs over the three p-like resonances listed above. We will use the form factor determined in 
Ref. jl8) in this work. The comparison with the experimental data of the ND E] and ARGUS collaborations J5l[] is 
shown in Fig. ^. 
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FIG. 9. The cross section for e + e 



The solid curve is based on the model of Ref. 



The experimental data are 



from the ND collaboration [J48| (circles) and the ARGUS collaboration |5lj (squares) 



C. pseudoscalar-axial vector meson annihilation 



In this class we consider chiefly irai — > 11, which is effectively a four pion process which we treat in the manner 
described previously. 

Thermal rate calculations indicate that this process is particularly important in the intermediate-mass region. 
There are, however, significant differences between calculations. Already in Ref. ^] it was shown that the dilepton 
production rates based on different models for the irpai dynamics can differ by an order of magnitude. This problem 
was recent revisited in Ref. ]5^] , where a comparative study was carried out for both on-shell properties and dilepton 
production rates using several models for the irpai dynamics. It was found that, although some of the models provide 
reasonable description of the on-shell properties, the corresponding dilepton rate could vary widely. By using some 
information from experimentally-constrained spectral function j53j, it was found that the effective chiral Lagrangian 
of Ref. [[34} , in which the vector mesons are introduced as massive Yang-Mills fields of the chiral symmetry, provide 
reasonable off-shell as well as on-shell properties of the irpai system. The interaction Lagrangian in this model is 
slightly elaborate and is given in Refs. [2q,K2[. The dilepton production cross section can then be obtained [B2[ 
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where TL is a complicated function of coupling constants, masses, and kinematics, and k v is the magnitude of the pion 
momentum in the center-of-mass frame. 
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FIG. 10. The cross section for e + e~ — > 7r + 7r~7r + 7r _ . The solid and dotted curves give our results in the first and second 
scenarios for the nai form factor, respectively. The experimental data are from the 772 collaboration (circles), [^Bj, the M3N 
collaboration [M (squares), and the ND collaboration |4Sj (triangles). 



One must now consider the issue of the electromagnetic form factor, IF^qJ. In principle, it can be determined 
by analyzing e + e~ — > 7r + 7r~7r + 7r~ and e + e~ — > 7r + 7r~7r°7r° data. Although many analysis have been carried out, 
an unambiguous determination of this form factor is not yet possible, since other intermediate-state can contribute 
to the same four-pion final state. Here we consider three scenarios, which should bracket this form factor. In the 
first scenario, we assume that most of the strength in the 7r + 7r~7r + 7r~ channel comes from a nai intermediate state. 
We can then determine \F 7Tai (M)\ 2 from experimental data for e + e~ — > 7r + 7r~7r + 7r~ from the 772 collaboration [j55| , 
the M3N collaboration |5(|, and the ND collaboration p§]. Further constraint on the IF^qJ 2 is provided by the 
experimental data for e + e~ — ► 7r + 7r _ 7r°7r°, which can also come from ir°uj intermediate state. This scenario basically 
sets an upper limit for the dilepton production from the na\ annihilation. In a second scenario, we assume that the 
7rai electromagnetic form factor is represented by the /?(770) only. We shall use the the parameterization of Eq. (4) for 
this purpose. The comparison of our results for e + e~ — > 7r + 7r~7r + 7r~ in the two scenarios with the experimental data 
is shown in Fig. |o|. Our cross section in the second scenario is similar to that extracted in Ref. |57|. In the second 
scenario, the missing strength is attributed to p(1700) which then directly decays into pixix without coupling to the 
■nai state. In Fig. we show the comparison of our results in the first scenario for the process e + e~ — > tt + tt^tt°tt° . 
The solid curve is the sum of our parametrizations for (7 e + e -^ 0l and cr e + e - ^-nu- Note that in this first scenario, the 
interpretations of the 47r ± and 27r ± 27r° data sets are consistent: our fit is adjusted to provide a compromise between 
the two sets of measurements. See Figs. 10 and |ll|. 
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FIG. 11. The cross section for e + e~ — > 7r + 7r - 7r°7r . The solid curve gives our fit in the first scenario for the ivai form factor. 
The experimental data are from the 772 collaboration, [^H), the M3N collaboration [^8| (squares), and the ND collaboration 
IH (triangles). 

Finally, as a lower bound we use the cross section cr e + e -_» Tai determined by the DM2 collaboration using partial 
wave analysis (PWA) J59|. We first parameterize the data as shown in Fig. [l^, and then use detailed-balance to obtain 
the cross section oVai-tiF- Note that the cross section at high masses is smaller than with our other two approaches. 
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FIG. 12. The cross section for e e — > na\. The open circles are the experimental data determined by the DM2 collaboration 
using partial wave analysis [p9| . The solid curve gives our fit to the data. 



III. THERMAL RATES 

The theoretical description of heavy-ion collision dynamics can be divided into two broad categories: transport 
(cascade) models and hydrodynamical models. In the latter model, (local) thermal and chemical equilibrium is 
usually assumed. To calculate dilepton spectra using the hydrodynamical model, one then needs to know the thermal 
dilepton emission rates, which we discuss in this section. 

According to kinetic theory, at a given temperature T, the differential rate per unit time and per unit volume for 
a meson pair to annihilate into a lepton pair, a + b — ► 11, is given by 

dR = dN r d 3 Pa d 3 Pb d 3 pi d 3 Pl 

dM 2 d A xdM 2 J 2E a (2n) 2 2E b (2Tr) 2 2Ei(2tt) 2 2Ej(2tt) 2 

x f(E a )f(E b )(2nf\M\ 2 S\ Pa + Pb - Pl - Pl )S(M 2 - (p, + Pl ) 2 ), (15) 

where / is the Bose-Einstein distributions for mesons, and M is an overall spin-isospin degeneracy factor. 
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FIG. 13. Dilepton emission rates from -kit and KK annihilation at three different temperatures. 



Approximating the Bose-Einstein distribution functions by the Boltzmann ones, we arrive at a simple expression 
for the thermal rate, 



Jtt rp 



(16) 



where Xi is a modified Bessel function, a(a + b — » ZZ) is the isospin-averaged cross section, and A is the usual kinematic 
triangle relation. 
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FIG. 14. Same as Fig. |l2j, for up and KK* + c.c. annihilation. 

We consider here three temperatures, T — 120, 150, and 180 MeV. The results for irir and KK are shown in Fig. 
|i~3| , the results for 7rp and KK* + c.c. are shown in Fig. |l4|, and the results for iroj and -nay are shown in Fig. |l5|. The 
results for irai are obtained in the scenario where the measured Air final states all come from an intermediate 7tai 
state. In the case of a\ meson, here we have also included the finite- width effects as in Ref. f60| . We mention that in 
the transport model, those are included for all the baryon and meson resonances. 
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FIG. 15. Same as Fig. |l^, for ncu and nai annihilation. 



The dilepton emission rates are seen to depend strongly on the temperature of the system. In heavy-ion collisions, 
as the system expands, its temperature rapidly decreases, so does the contribution to the dilepton rate. So most of 
the dileptons are emitted in the early hot and dense stage of the collisions. Comparing with previous calculations, 
our results for ttit and KK are in agreement with those of Ref. [p7| , and our results for -np and KK* + c.c are in 
agreement with those of Ref. |2^| . Finally, as expected, our results for ttlu and nai are different from those in Refs. 
pq , |60[ , since we use new form factors for these processes which we obtain from e + e~ annihilation data. Of all the 
processes considered here, the nai is found to be the most important one for intermediate-mass dileptons. This is in 
line with previous observations p8| , ^0|j60| . 

Dilepton emission rates from hot and dense hadronic matter have been calculated in different approaches. Here we 
discuss and compare our results with three of these approaches: the spectral function approach of Ref. [ |53| |, the chiral 
reduction formalism of Refs. |6l|]6lf |, and a model based on quark- hadron duality [|S3| . 

In general, the dilepton emission rate is given by the thermal expectation value of the electromagnetic current- 
current correlation function |^, |6^ , |65|| , 



dR 

dM 2 



6tt j M 2 



d*xcr i9X ((ji t (x)j M (o)))i 



(17) 



Different approaches differ in the way this correlation function is approximated and calculated. In Ref. |53|| , this 
current-current correlator was reduced to a number of spectral functions which can be extracted from the e + e~ 
annihilation and r lepton decay data. Including the soft final-state interaction corrections that leads to the parity 
mixing phenomenon (i.e., the mixing between the vector and the axial- vector currents), the dilepton emission rate 
can be expressed as [B3| 
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<II{ MTK^M/T) (p* m {M) - (e- j)(p v (M) - /(M))) , (18) 



dM 2 



where e = T 2 /6F 2 , and p em , p v , and p" 4 are the electromagnetic spectral function, the vector spectral function, 
and the axial-vector spectral function, respectively. p em can be related to the total hadronic cross section of e + e~ 
annihilation, 

p e m{M) = M 2 a(e + e-^hadrons) 

167T J Q! 2 

Similarly, the vector spectral function can be determined by selecting events with even number of pions 

M 2 

p v ( M ) = E ^ e+e " -> 2mr )- ( 2 °) 

n— 1 

Finally, the axial- vector spectra function p" 4 can be extracted from the differential probability of the r lepton decaying 
into an odd number of pions. Thus, with these spectral functions extracted from the e + e~ and r decay data, one can 
determine the empirical dilepton emission rate using Eq. (O). Note that the terms in e represent a small correction. 

In Ref. HI , the idea of a Hagedorn resonance gas approach was explored. There, the vector meson mass distribution 
is characterized by the mass spectrum py(m) (we use p v for the vector spectral function, and pv for the vector 
meson mass spectrum). The dileptons are associated with direct vector meson decays (the vector meson dominance 
assumption). Therefore the thermal rate can expressed as 

Pv{M)- — < ^——M 2 TK 1 (M/T), (21) 



dM 2 rvy Qitg 2 (M) 

where g(M) is the vector- meson-photon coupling constant, which also determines the vector meson production cross 
section in the e + e~ annihilation, 

In Ref. Q, it was assumed that the total cross section of e + e~ annihilation is saturated by the production of vector 
mesons, which is supported by experimental observation p3| , 

a(e+e- -> hadron) = R^pJL — = CT ( e + e " _> y), (23) 

where R ex P is the experimental R value, which reflects the strong interaction aspect of the e + e _ annihilation physics. 
With this approximation the dilepton emission rate can then be written as 

^=R^MT Kl (M/T). (24) 

This is very similar to the results of Ref. ]53|], except here that the soft-pion corrections are neglected. In Ref. |33|, it 
was further assumed that R ex P can be approximated by RP art = N c J2f e /j the parton model prediction for R. The 
dilepton emission rate is then determined theoretically without referring to experimental data, 

J^ =RP ar t ^_ MTKl(M/n (25) 

For three flavor case, RP art = 3(4/9 + 2/9 + 2/9) = 8/3. This equality states that the dilepton emission rate from 
hadronic gas can be approximated by the parton model prediction with the same accuracy as the parton model 
describes the total hadronic production cross section in e + e~ annihilation. This approximation would apply to 
invariant masses greater than ss 1.5 GeV. 

Finally, in Refs. |6l| , ^2| , a chiral reduction formalism f36|| was use to reduce the current-current correlation function 
in Eq. ( |l7| ) into a number of vacuum correlation functions, which are constrained by the experimental data from 
e + e _ annihilation, r lepton decay, pion radiative decay, and two-photon fusion reactions. This approach allows for a 
systematic expansion in terms of pion density (and nucleon density if baryons are included) . 
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FIG. 16. Dilepton emission rates from different approaches. 



In Fig. 16, we compare the dilepton emission rates obtained with our different scenarios with those obtained in the 
three models just outlined. The solid, long-dashed, and dash-dotted curves are our results based on our pictures for 
the 7rai cross sections. The scenario that includes the p(1700) in the irai form factor leads to a larger rate above 1.2 
GeV invariant mass, as expected. The open circles are the results of Ref. |33| based on empirical spectral functions. 
Our result with a full 7rai form factor is in agreement with the spectral function result up to 1.6 GeV invariant mass. 
It is slightly above at higher invariant masses. The scenario that excludes p(1700) in the nai form factor is slightly 
lower than the spectral function result between 1.2 and 1.6 GeV invariant mass. In other mass regions, they give 
comparable rates. Finally, the use of the DM2 partial wave analysis for the irai cross section yields rates that are 
somewhat lower than those obtained in the other models, passed 1.75 GeV invariant mass. We will see later that hard 
background processes will mask this difference. 

The short-dashed curve gives the results of Ref. which is based on a SU(3) version of the chiral reduction 
formalism. The result is in good agreement with the one from empirical spectral functions |53| ], suggesting that 
higher-order interactions (not explicitly included in Ref. |j53f ) are not significant at the temperatures considered here. 
Finally, the dotted curve gives the result of Eq. ( p5| ) which makes use of the hadron-quark duality that allows the 
use of RP art in the place of R exp . It is interesting to note that the dilepton rate from this approach is in very good 
agreement with the one obtained with empirical spectral functions above about 1 GeV. This supports the idea that, 
even though there are resonances in the the intermediate mass region, these resonances cannot be easily resolved 
individually. 

As a partial summary, our results are based on kinetic theory which assumes pair-wise collisions of mesons. To 
make use of e + e~ annihilation data as a constraint, we need to know the hadronic production cross section in specific 
channels, such as tttt, KK, or nai. Data are available for most of the processes considered here. But for nai an 
unambiguous determination of its production cross section from the measured An cross sections is not yet possible. 
This introduces model dependence. On the other hand, the other three approaches use directly the total hadronic 
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production cross section in e + e~ annihilation, without the need to separate them into specific channels and thus 
avoiding the uncertainties involved in such a separation. The nice agreement shown in Fig. [l6] thus suggests that 
the model dependence in our approach is not severe. Furthermore, this level of agreement from different approaches 
for these elementary processes provides us with the basis for further discussions about the dynamics of heavy-ion 
collisions, the possibility of GQP formation, and charm enhancement. 

IV. APPLICATION TO HEAVY-ION COLLISIONS AT SPS ENERGIES 

To compare with experimental data, the thermal rate calculation is certainly not sufficient. One needs a transport 
model that describes the dynamical evolution of the colliding system, and integrate the dilepton production rate over 
the entire reaction volume and time. In heavy-ion collisions at CERN SPS energies, many hadrons are produced 
in the initial nucleon-nucleon interactions. This is usually modeled by the fragmentation of strings, which are the 
chromoelectric flux-tubes excited from the interacting quarks. One empirically successful model for taking into 
account this nonequilibrium dynamics is RQMD |j67f . To define a relativistic transport model for heavy- ion collisions 
at these energies, we have used as initial conditions the hadron abundances and distributions obtained from string 
fragmentation. Further interactions and decays of these hadrons are then taken into account as in most relativistic 
transport models. This approach is found to provide a good description of hadronic observables (such as transverse 
mass spectra and rapidity distributions) in heavy-ion collisions at CERN SPS energies fri||26|| . 

We should make here some comments on the treatment of resonances and multi-particle interactions in the transport 
model. The narrow-width approximation was not used in our calculation. The mass distribution of the resonances 
(for example the a\) in our approach is taken into account. Note that this realistic sampling of the spectral function 
is different from usual kinetic theory estimates of the dilepton emission rates where the resonance is usually given a 
fixed mass. Furthermore, in our transport model, the fact that the resonance has a finite life time is fully respected: 
its population is evolved dynamically. The resonance is formed in collisions and also eventually decays. In the early 
stage of the heavy-ion collision, resonance formation is favored since the density is high. As the system expands and 
the density gets smaller, the resonance decay process gradually becomes dominant. This is again different from kinetic 
theory calculations. 

Let's take it + a\ — > dilepton as a concrete example. In our approach, a rho meson is first made from two-pions 
colliding. It then interacts with a third pion to form an a\ meson. Finally this a\ meson annihilates a fourth pion to 
produce a lepton pair. The explicit treatment of intermediate resonances takes care of the formation time typical of 
strong interaction, which is very important for heavy-ion collisions at SPS energies. Furthermore, this algorithm also 
includes the possibility that in the course of the multi-step process the intermediate rho or a\ meson might interact 
with other particles, say a nucleon, rather than with the third or the fourth pion. This amounts to a dynamical 
generation of a broadened width. Finally, for any reasonable comparison with experimental data from heavy-ion 
collisions, a transport model that describes the dynamic evolution of the colliding system is a required tool. 

The contributions from the secondary processes outlined above are shown in Fig. |l7]. These are obtained in our 
relativistic transport model |l7|j26| ], including the HELIOS-3 acceptances, mass resolution, and normalization JlO| . 
Here N^^/Nch represents number of dimuon pairs within the 50 MeV mass bin, normalized by the charged particle 
multiplicity. It is seen that the nai process is by far the most important source for dimuon yields in this mass region. 
The ttu process also plays some role in the entire intermediate-mass region, while the contributions from mr, irp and 
KK are important around 1 GeV invariant mass The relative importance of various processes follow basically the 
thermal rate predictions, as discussed in the previous section. 
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FIG. 17. Dimuon spectra in central S+W collisions at 200A GeV from different secondary processes. 



In Fig. [18|, we add the secondary contributions obtained in our transport model to the background, and compare 
again with the HELIOS-3 data for central S+W collisions. It is seen that the data can be adequately reproduced. 
This highlights for the first time the importance of the secondary processes for the intermediate-mass dilepton spectra 
in heavy-ion collisions. This is an important step forward in the use of intermediate-mass dilepton spectra as a 
probe of the phase transition and QGP formation. Although the current data do not show any necessity to invoke 
the QGP formation in S-induced reactions, consistent with some conclusions from J/'F physics, we believe that 
the observation that the secondary processes do play a significant role in the intermediate-mass dilepton spectra is 
interesting and important. We note in passing that the change of slope observed in the experimental data corresponds 
in our interpretation to a crossover between the secondary processes and the (Drell-Yan/charm) background. This 
represents a transition from soft to hard physics and this observation should be explored further in the future. 
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FIG. 18. Dimuon spectra in central S+W collisions at 200A GeV. The dotted curve gives the background contribution as 
shown in Fig. [j], the dashed curve gives the contribution from secondary processes shown in Fig. [j^, and the solid curve gives 
the sum of the two contributions. 



The results shown in Figs. |l7|and[l§| are obtained under the assumption that all the 4tt final state in the e + e~ cross 
section proceeds through the 7tai intermediate state. We also did a calculation in which the 7tai form factor contains 
only the normal p(770). The results are shown in Fig. [H| by the dotted curve. Finally, to complete our survey of 
possible constraints and uncertainties in nai — > e + e~ cross sections, we did a calculation in which this cross section 
is obtained from the e + e~ — > nai cross section determined by the DM2 collaboration in partial wave analysis (PWA) 
p9[ (see Fig. |l2]). The results are shown by the dashed curve in Fig. |l9|. The region between the solid and dashed 
curves thus reflect the uncertainty for heavy-ion collisions due to our limited knowledge of the 7rai cross section. This 
area is not unreasonably large. From a formal point of view, it is fair to say that no strong evidence currently exists 
coupling the p(1700) to a 7rai state f68|| , even though better Air data could help resolve this issue along with others of 
interference and irhi contribution. Note that the fact that the DM2 results are lower than in other scenarios at high 
invariant masses is concealed by the hard background. 

Also shown in this figure is the results of Ref. [ ]62| which uses the dilepton emission rate obtained in the chiral 
reduction formalism with a simple fireball model for the dynamics of heavy-ion collisions. It is see that their results 
are also in good agreement with the experimental data. 
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FIG. 19. Dimuon spectra in central S+W collisions at 200A GeV. The solid and dotted curves give our results with and 
without p(1700) in the na\ form factor, while the dashed curve gives our results using DM2 PWA cross section. The dashed 
curve gives the results from J32|. 



Another issue we want to address here is the effects of dropping vector meson masses on the entire dimuon spectra 
from the threshold to about 2.5 GeV. In Ref. |l7j it was shown that the enhancement of low-mass dileptons could be 
interpreted as a signature of vector meson masses decreasing with increasing density and temperature. This should 
affect the dilepton spectra in the intermediate-mass region, mainly through two effects. One is the change of the 
invariant energy spectra of these secondary meson pairs. The second effect enters through the modification of the 
electromagnetic form factor. Since we can only conjecture how the masses of the higher vector resonances change 
with density and temperature, we shall assume for simplicity that they experience the same amount of scalar field as 
the "common" rho meson, namely, my V , = myy — 2/3g a {a) ]l7| j. The results of this calculation are shown in Fig. 
p0| . Below 1.1 GeV and especially from 0.4 to 0.6 GeV, the agreement with the experimental data is much better 
when the dropping vector meson mass scenario is introduced, as was already shown in Ref. [ pT| . At higher masses, the 
dropping mass scenarios somewhat underestimates the experimental data. Also note that our starting point was our 
first scenario. For completeness, however, we have to state that there might be additional contributions from, e.g., 
secondary Drell-Yan processes |69(| that were not included in this study. There could also be enhanced production 
of charmed mesons. Furthermore, as we progress higher in invariant mass, the role of baryons has to be carefully 
assessed. So far, the baryons seem to play little role in the overall dilepton yield [|7(J. This statement was examined 
in Re f. |6l|] for masses below 1 GeV, which is being extended to the intermediate-mass region in |p2| . See however 
Ref. for an alternate discussion on the role of baryons. 
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V. SUMMARY AND OUTLOOK 



In summary, we have carried out a study of intermediate-mass (between 1 and 2.5 GeV) dilepton spectra from 
hadronic interactions in heavy-ion collisions. The processes included are ttt: — ► 11, up — » II, nai — ► ZZ, 7ro> — » ZZ, 
ifX — » ZZ~, and KK* + c.c — > ZZ. We calculated the elementary cross sections for these processes based on chiral 
hadronic Lagrangians for pseudoscalar, vector, and axial-vector mesons. The corresponding electromagnetic form 
factors are determined by fitting to the experimental data for the reverse processes of e + e~ — > hadrons. While 
the form factors for irir, KK, irp, iru>, and KK* + c.c. are uniquely constrained by the e + e~ data, there are still 
uncertainties in the ttqi form factor, for which we have considered three scenarios. 
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FIG. 20. Dimuon spectra in central S+W collisions at 200A GeV. The solid and dotted curves give our results with and 
without meson medium effects. 



We used these cross sections with kinetic theory to calculate the dilepton emission rates at finite temperature. The 
■na\ — ► 11 process was found to be the most important one in the intermediate-mass region, as was found in previous 
calculations. Our rates were compared with those from other approaches. This provides us with the confidence that 
we have some control of the dilepton spectra from the secondaries hadronic interaction. This is important for the use 
of the dilepton spectra as a probe of phase transition and QGP properties. 

Finally we apply these elementary cross sections in a relativistic transport model and calculate dilepton spectra in 
S+W collisions at SPS energies. Again, the nai — ► 11 was found to be the most important process. The comparison 
of our results with experimental data from the HELIOS-3 collaboration indicates the importance of the secondary 
hadronic contributions to the intermediate-mass dilepton spectra. 

This work can be extended in several directions. Although the current experimental data from the HELIOS-3 
collaboration arc not very sensitive to the different assumptions of the ixa\ form factor, it would be very useful to 
pin this down since the nai is the single most important source of the intermediate-mass dileptons from hadronic 
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secondary collisions. Also, a more elaborate calculation of the charmed meson production and initial (as well as possible 
secondary) Drell-Yan contributions in heavy-ion collisions is needed. Such a calculation could put the conclusions 
reached in this work concerning the role of secondary hadronic processes and dropping vector meson masses on an even 
firmer footing. Finally, the current investigation can be extended to higher incident energies, such as those of the RHIC 
collider, by combining the cross sections (or thermal rates) obtained in this study with, e.g., hydrodynamical models 
for the evolution of heavy-ion collisions at the RHIC energies. This kind of study is useful for the determination 
of hadronic sources in the dilepton spectra to be measured by the PHENIX collaboration, and thus for the clear 
identification of the dilepton yield arising from the QGP. 

We are grateful to Nikolai Achasov, Gerry Brown, Sandy Donnachie, Axel Drees, Madappa Prakash, Ralf Rapp, 
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